Symmetry Breaking on the Three-Dimensional Hyperkagome Lattice of Na4lr308 
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We study the antiferromagnetic spin- 1/2 Heisenberg model on the highly frustrated, three-dimensional, hy- 
perkagome lattice of Na4lr308 using a series expansion method. We propose a valence bond crystal with a 72 
site unit cell as a ground state that supports many, very low lying, singlet excitations. Low energy spinons and 
triplons are confined to emergent lower-dimensional motifs. Here, and for analogous kagome and pyrochlore 
states, we suggest finite temperature signatures, including an Ising transition, in the magnetic specific heat due 
to a multistep breaking of discrete symmetries. 
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Introduction. — In the last decades, the search for exotic 
magnetic states has mainly concentrated on low-dimensional 
systems, d < 2, where the effects of fluctuations are partic- 
ularly strong [1-4]. It has only been realised more recently 
that — whereas magnetic order is more easily destabilised in 
low dimension — certain unusual phases are in fact more ro- 
bust in d = 3 [5]. 

A case currently under active investigation is the compound 
Na4lr308 [6], where a subtle interplay of spin and orbital de- 
grees of freedom leads to an effective S = 1/2 description 
[7] of the magnetic Ir ions occupying the sites of a hyperk- 
agome lattice of corner- sharing triangles (Fig. 1). This ma- 
terial shows no sign of magnetic ordering down to the lowest 
temperatures, thus presenting a candidate for a quantum spin 
liquid [8-10]. Here, we follow a different lead, namely the 
possibility of the magnetic state being a valence bond crystal 
(VBC), breaking translational but not spin-rotational invari- 
ance. Interestingly the gadolinium gallium garnet lattice [11] 
with nearest neighbor interactions consists of two indepen- 
dent hyperkagome networks of opposite chirality, providing a 
broader scope for our VBC proposal. 

We find a candidate ground state with a huge (in this case 
72-site) unit cell. We elucidate its structure and discuss its 
excitations, which include emergent low-dimensional motifs. 
We discuss the general implications of these large-unit-cell 
VBCs for the expected (sequence of) phase transitions as 
the temperature is lowered. These considerations also ap- 
ply to analogous states proposed for kagome [12-14] and py- 
rochlore [15-18] magnets. We also address some general fea- 
tures of the series expansion method [12, 13] used to derive 
our effective Hamiltonian. 

Lattice structure and effective model — The hyperk- 
agome lattice [Fig. 1(a)] is a three-dimensional lattice of cor- 
ner sharing triangles and as such a cousin of the kagome lat- 
tice in d = 2 shown in Fig. 1(b). It can be obtained from 
the pyrochlore lattice by removing one site per tetrahedron, 
yielding a connected network of triangles with a unit cell of 
twelve sites in which all sites are equivalent [6]. We set the 
length of the cubic unit cell to one in the following. An im- 
portant property is that the shortest cycle (closed loop) on this 
lattice beyond the triangles involves ten bonds [cf. Fig. 1(c)] 
and will be called "decagon" in the following, in analogy with 




Figure 1 : (Color online) Lattice structure: (a) cubic unit cell of the 
hyperkagome lattice (thin bonds and small sites show the extra bonds 
and sites of the pyrochlore lattice), (b) the 2D kagome lattice, (c) a 
'decagon', the smallest loop of the hyperkagome lattice, highlighting 
the structure of two intersecting planes by different coloring. 



the length six "hexagon" loops on kagome. 

We aim to determine the quantum mechanical ground state 
of the antiferromagnetic nearest-neighbor S = 1/2 Heisen- 
berg model on the hyperkagome lattice with N sites and peri- 
odic boundary conditions given by the Hamiltonian: 

H = JY,S,-Sj . (1) 

Focusing on valence bond crystal formation, we perform 
series expansions around arbitrary nearest-neighbor singlet 
(dimer) coverings of the hyperkagome lattice to obtain an ef- 
fective Hamiltonian (energy function) for the 2^/^+^ dimer 
coverings, following earlier work of Singh and Huse on the 
kagome lattice [12, 13]. We subsequently determine and char- 
acterize the dimer coverings minimizing the energy. While the 
triangles covered by a dimer are local ground states of Hamil- 
tonian (1), the N/6 triangles not hosting a dimer induce quan- 
tum fluctuations and if strong enough can melt the dimerised 
state. 

We start with an arbitrary dimer covering, V, of 'strong' 
singlet bonds and treat the remaining part of as a perturba- 
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This procedure folds the quantum fluctuations of the empty 
triangles into effective potential energy terms (given by power 
series in A) for various local motifs appearing in the dimer 
coverings. If a VBC is actually favored — or is at least locally 
stable — such an expansion should be well behaved, converg- 
ing towards the correct ground state and energy thereof at rea- 
sonably low expansion orders also at A = 1. 

In our analysis, we expand the Hamiltonian up to fifth or- 
der in A. To this order, the Hamiltonian consists of a con- 
stant term, a term favoring two empty triangles linked by one 
dimer and a term favoring five empty triangles arranged along 
decagons. The effective model is illustrated in the following: 
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(2) 

The coefficients are given as a = 9AV32 + 15A'^/2048 + 
AV128, (3 = 5A'^/264 + AV2048 and 7 = 105AV8192. 
Note that since the number of empty triangles is the same for 
all dimer coverings, the term proportional to a simply corre- 
sponds to a global energy shift. 

Minimization of the effective Hamiltonian. — Given the 
constrained nature of dimer coverings, minimizing the clas- 
sical dimer Hamiltonian (2) among all Nx) = 2^/^+^ dimer 
coverings of the hyperkagome lattice is a non-trivial task. 
In order to find ground states we resort to a classical dimer 
Monte Carlo approach based on non-local worm updates [19]. 
Simulating the effective classical dimer model (2) at A = 1 
on a range of cubic samples of = x 12 sites with 
L = 3, 4, . . . , 9, 12, 15, 18 we found states which reach the 
lower bound for the energy of the term binding two neigh- 
boring empty triangles (an empty triangle can have at most 
three neighboring empty triangles). On an L = 6 sample we 
uncovered our lowest energy state state (to be described be- 
low) which satisfies all two-triangle interactions and has 144 
satisfied decagon terms. Simulations involving only decagon 
energy terms reached a slightly larger number of 147 perfect 
decagons on the same sample. So it seems — unlike the case 
of the kagome lattice — that there is no way of minimizing the 
shortest closed loop term and the neighboring empty triangle 
interaction simultaneously on the hyperkagome lattice. Given 
that the energy penalty is quite small (3 out of 147) and the 
decagon coupling 7 is smaller than the empty triangle binding 
term /3, we are optimistic that the reported VBC state is the 
bona fide ground state of the classical dimer model Eq. (2). 

As shown in Table I the candidate VBC on the hyperk- 
agome does very well in minimizing the energy as compared 
to the best spin liquid estimate given in the literature [10]. 
In fact, it seems to compare even more favorably to the spin 
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Table I: VBC energies per site in units of J for the infinite system 
evaluated at A = 1 for the hyperkagome versus the lowest spin liquid 
(SL) energy quoted in the literature [10]. 



liquid than is the case in two dimensions (comparing the re- 
sults obtained in Refs. 13, 20). On general grounds one may 
also expect that a VBC state is likely to be more stable on a 
lattice with longer closed loops. To higher order in A, increas- 
ing numbers of terms involving geometrical motifs growing 
in size will enter the effective Hamiltonian, which will generi- 
cally compete with each other. We note that sub-leading terms 
tend to increase the magnitude of the coefficient of pairwise 
empty triangle interactions, (3, while the coefficient of the 
closed loop terms, 7, decreases at sub-leading order. Thus, 
it is likely that the regime where our candidate state wins is 
maintained at higher orders. 

We have also established that the series expansion itself has 
the capability of signaling its own breakdown by explicitly 
carrying out the analysis on the checkerboard lattice where 
long range quadrumer order (rather than dimer order) is ex- 
pected [21]. In this case the series favor different ground state 
symmetries already at the two first non-trivial orders (second 
and third in this case) when evaluated at A = 1. This is in 
sharp contrast to the kagome and hyperkagome expansions. 
Moreover, it was demonstrated explicitly that the energy of 
the candidate VBC changes only little even between fifth [12] 
and eighth [13] order on the kagome lattice. 

Description of the VBC state. — Our proposed VBC state 
has a unit cell of 72 sites [30] (containing 36 dimers) and is 
illustrated in Eig. 2(a). The two main building blocks are: 
i) two kinds of decagon-rods running along a [110] direction 
each consisting of a chain of stacked perfect decagons, with 
two neighboring decagons sharing a dimer and two empty tri- 
angles [cf. Eig. 2(c)] and ii) A-chains — with two degenerate 
dimer coverings for each chain — running parallel to the rods, 
see Eig. 2(d). There is an equal number of decagon rods and 
A-chains in the VBC structure, as can be inferred from their 
mutual arrangement in the [110] planar projection displayed 
in Eig. 2(b). The two building blocks are reminiscent of i) the 
perfect hexagon network and ii) the pin wheels of the Marston- 
Zeng VBC [12, 14] on the kagome lattice. 

The 36 dimers per unit cell of the crystal can be assigned as 
follows: 18 form links between empty triangles, out of which 
16 belong to the two decagon-rods while two connect rods in 
only one of the two directions transverse to the rods, forming 
two dimensional layers within in the three dimensional lattice 
as shown in Eigs. 2(a,b). Out of the remaining 18 dimers, 
12 are located on the A-chains, while the other 6 shield the 
A-chains from the network of empty triangles located on the 
decagon rods. 

A decagon-rod can be formed along a backbone [indicated 
by the line of blue and red bonds in Eig. 2(c)] in either of 
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Figure 2: (Color online) VBC structure: (a) details of the VBC ground state on a slab of the hyperkagome viewed along a [110] direction: 
dimers are red and weak bonds blue except on the decagon-rods, which are entirely yellow. The empty triangles are indicated by green 
spheres. The two-dimensional, kagome-like, projection of the state is shown in (b), highlighting the positions and orientations of the decagon- 
rods (yellow) and the A-chains (red vs green). The shaded region corresponds to one choice of a VBC unit cell. The positions of the backbones 
of the decagon-rods are marked by black and their bridging dimers are indicated in brown, (c) shows a decagon rod — the main building block 
of the VBC ground state. The decagon-rods are directed along a backbone (blue). The two degenerate dimer coverings of a A-chain (red vs 
green) are shown in (d). 



the six equivalent [110] directions inherited from the straight 
lines of bonds of the underlying pyrochlore lattice before de- 
pletion. It turns out that choosing a backbone and one of two 
possible dimer covering thereof fixes the rod and subsequently 
the entire VBC completely (up to the individual two-fold de- 
generacy of the A-chains), which implies a total of 72 distinct 
VBC ground states. 

Finite temperature ordering transitions — The VBC we 
have described is of course not the first one to be proposed for 
frustrated quantum magnets which exhibits a large unit cell 
with a complex ordering pattern — see proposals for kagome 
[14, 22] or pyrochlore [15-18]. How is one to verify or falsify 
these proposals? This is all the more pressing since one is 
unlikely to encounter the 'perfect' Hamiltonian, Eq. (1); nor 
will it be easy to probe any given material below a potentially 
tiny ordering temperature. Can one nonetheless make some 
qualitative predictions? 

With this motivation, let us discuss the behaviour of the 
magnet as T is raised. Whereas it is in principle conceivable 
for there to be a single transition into the paramagnet, this 
is rather unlikely: there is no a priori reason why one of the 
different components of the complex ordered state should not 
melt before the others, leaving a partially ordered state be- 
hind. Our VBC involves for instance a Z2 degree of freedom 
encoding which sublattice of triangles hosts the triangles not 
occupied by a dimer (the triangles of the hyperkagome lattice 
exhibiting a bipartite structure). The expansion in A in fact 
suggests that this is the most robust feature, predicting a rel- 
atively high temperature Ising transition separated from other 
transition(s) which occur at lower temperature [31]. The situ- 
ation in the kagome [14] and pyrochlore [15-18] examples is 
entirely analogous. 

Excitations and emergent low -dimensional structures. — 
The proposed VBC state features a huge number of low-lying 
singlet excitations. For instance, we expect singlet excitations 



involving local rearrangements of dimer motifs living at an 
energy scale of the order of the couplings /3J^1.9xlO~^J 
and 7J ^ 1.3 X 10~^ J at A = 1. A subextensive number 
corresponding to the two distinct arrangements of dimers on 
each A-chain remains degenerate until even higher order in 
perturbation theory. 

The natural triplet excitation picture at small A is to excite a 
strong dimer from the singlet to a triplet state (called triplon) 
with an energy cost J. In first order perturbation theory this 
triplon will possibly delocalize by hopping processes between 
strong dimers. Analyzing the resulting tight-binding model of 
triplons we find that triplons on the A-chains and the shielding 
dimers (18 dimers in total) remain completely localized on 
their respective dimer due to the absence of hopping terms in 
first order. The remaining 18 dimers in the VBC unit cell form 
12 dispersive and 6 (lowest lying) flat triplon bands confined 
to a two-dimensional layer of decagon rods and their bridging 
bonds. The resulting dispersion is shown in Fig. 3. Of the six 
triplons in the flat band (shown in brown in Fig. 3), four may 
be localized on the four perfect decagons in the VBC unit cell, 
while the remaining two on loops with 8 dimers involving the 
bridging dimers and extending beyond a single unit cell. 

When going beyond the first order treatment we anticipate 
similar irregular behavior of the series as on the kagome lat- 
tice [13]. At any rate, the huge number of low lying excita- 
tions makes it seem very likely that coherent excitations will 
exist only in a very small window in energy-momentum space. 

A complementary magnetic excitation picture has recently 
been developed in Ref. 23, based on the presence of gapped, 
deconfined spinons on the Husimi cactus (Bethe-lattice ana- 
logue of the kagome and hyperkagome lattice) and the one- 
dimensional A-chain. Since our VBC state explicitly features 
rather isolated A-chains it is natural to assume that the com- 
pletely localized triplons obtained in the first order treatment 
above actually decay into a spinon continuum when approach- 
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Figure 3: (Color online) Dispersion of the triplons to first order in 
A, plotted as function of the momenta ka and kb in the decagon- 
rod plane spanned by the vectors a and b defined in Fig. 2. The 
dispersion is flat in the perpendicular kc direction for all the modes. 
The inert triplets at constant energy J are not shown. 

ing the homogeneous A = 1 limit. 

Numerical simulations of the dynamical spin response on 
the kagome lattice [24] suggest that a combination of both 
pictures might be appropriate. At energies just above the spin 
gap the spectral density is quite high, reminiscent of the coher- 
ent response of triplon modes, while at higher energies a rather 
featureless continuum is visible, reminiscent of the deconfined 
spinon picture. The confinement of triplon excitations to two- 
dimensional decagon-rod planes suggest that spectral function 
features are almost momentum independent in the direction 
transverse to these planes. 

Discussion and relevance to experiment. — In this Letter 
we have presented a candidate valence bond crystal ground 
state for the spin- 1/2 Heisenberg model on the hyperkagome 
lattice. The regularity and favorable energetics of our series 
expansion establishes the VBC as a serious contender to the 
earlier spin liquid proposals [8-10]. It would be interesting 
to investigate this further using the "resonating" approach of 
Ref. 25, which treats quantum dynamics in a different manner. 

To what extent a Heisenberg model is relevant for low en- 
ergy physics of the recently studied Na4lr308 material is a 
matter of current research — its relevance has been argued for 
by several authors (see eg Refs. 6, 7, 26) and assumed by oth- 
ers, while there are also reasons to believe that the material is a 
weak Mott insulator undergoing a transition to a metal under 
pressure, thus making the assumption of a nearest-neighbor 
Heisenberg description thereof doubtful (see eg Refs. 27- 
29). However, we note that the lack of magnetic ordering 
as well as the existence of a peak in Cy/T in the exper- 
iment is a priori consistent with our picture. We also sug- 
gest that an anisotropic spin structure factor arising from the 



low-dimensional excitations discussed above serve as a fur- 
ther probe of our theory. 

We acknowledge useful discussions with F. Alet, Z. Hao, 
D. Huse, K. Penc, I. Rousochatzakis and H. Takagi. 
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